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In the quest to miniaturise photonics, it is of paramount importance to control Ught at the nanoscale. We 
reveal the main physical mechanism responsible for operation of gap plasmon-based gradient metasurfaces, 
comprising a periodic arrangement of metal nanobricks, and suggest that two degrees of freedom in the 
nanobrick geometry allow one to independently control the reflection phases of orthogonal light 
polarisations. We demonstrate, both theoretically and experimentally, how orthogonal linear polarisations 
of light at wavelengths close to 800 nm can be manipulated independently, efficiently and in a broad 
wavelength range by realising polarisation beam splitters and polarisation-independent beam steering, 
showing at the same time the robustness of metasurface designs towards fabrication tolerances. The 
presented approach establishes a new class of compact optical components, viz., plasmonic metasurfaces 
with controlled gradient birefringence, with no dielectric counterparts. It can straightforwardly be adapted 
to realise new optical components with hitherto inaccessible functionalities. 
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In the ongoing quest to miniaturise photonics to the micro- and nanoscale, conventional optical components, 
such as dielectric lenses and polarisers, are no longer adequate for steering and manipulating light owing to 
their inherently bulky nature. One of currently investigated approaches to resolve this issue is the application 
of metamaterials, which allow one to design materials with specific and exotic optical properties, including 
negative refraction^'^ and moulding the flow of light using inhomogeneous metamaterials and transformation 
optics^'^. 

An alternative approach towards light control at the nanoscale, which also seems more feasible for practical 
realisation and implementation, is concerned with two-dimensional (2D) metal nanostructures with sub wave- 
length periodicity, also known as metasurfaces. Despite having negligible thicknesses as compared to the incident 
wavelength, metasurfaces are characterised by the ability to strongly manipulate both the amplitude and phase of 
impinging light near (plasmonic) resonances of the unit cell constituents^. It should be however noted that 
resonant elements allow one to control the phase only in a limited range, from 0 to tt (radians), due to their 
Lorentz-like polarisabilities^"^. In order to access the whole phase space of 2n, two different approaches are 
typically employed. In the first approach, the 27i-space is achieved with transmitted/reflected light polarised 
orthogonally to the incident wave by using V-shaped nanoantennas that can be designed to realise a variety of flat 
optical components^""^^. Nevertheless, the cross -polarisation condition sets an upper bound to the efficiency of 
the devices that in the limit of negligible absorption is theoretically estimated to 25%^^. Importantly, all demon- 
strated functionalities can also be realised, in principle, by profiling the surface of a dielectric substrate, e.g., with 
ion-beam milling technique, in the fashion employed for blazed gratings^^. The second approach is similar, at least 
from the first sight, to that used in reflectarray antennas^^ in which an array of metal rectangular patches is placed 
in close proximity to a metal film or substrate, only separated by an optically thin dielectric spacer. The metal- 
backed design ensures an almost full control of the phase space together with a high efficiency, reaching 100% in 
the limit of negligible Ohmic losses. These two attractive properties have been demonstrated in the realisation of 
efficient wave plates^^"^\ focusing mirrors^^'^^, blazed gratings^^, and components interfacing propagating and 
surface waves^l In the aforementioned work, the metasurface functionalities are related to a proper design of the 
phase profile for one polarisation, mimicking the phase accumulation acquired when light transverses conven- 
tional optical components. Metal-backed metasurfaces, however, seem to be much more versatile according to 
very recent numerical simulations of metasurfaces showing polarisation- sensitive reflection in the near-infrared^^ 
and infrared^^ regimes. Note, though, that a design technique and experimental verification of such gradient 
birefringent metasurfaces are still lacking. 

Here, we elucidate the underlying physics responsible for operation of metal-backed metasurfaces in visible 
and near-infrared, emphasising the importance of excitation of gap surface plasmons (GSP) under metal patches 
or nanobricks^°'^^'^l The developed physical understanding allows us to generalise the design rules of plasmonic, 
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Figure 1 | Independent control of TE and TM polarisation, (a) Sketch of basic unit cell consisting of a gold nanobrick on top of a glass spacer and gold 
substrate, (b) Calculated reflection coefficient ras a function of nanobrick widths for A = 240 nm, t = ^5 = 50 nm, and X = 800 nm (for further details see 
Supporting Information). Colour map shows the reflection coefficient amplitude for TM polarisation, while lines are contours of the reflection phase for 
both TM and TE polarisation, (c) Super cell of metasurface functioning as a polarisation beam splitter (case 1 ): TM (TE) waves are reflected into +1 ( — 1 ) 
diffraction order. (d,e) Theoretical performance of the metasurface for TM polarisation, displaying (d) the x-component of reflected E-field (EJJ) 
just above the metasurface at the design wavelength (the amplitude of incident E-field is 1 V/m), and (e) amount of incident light reflected into the lowest 
diffraction orders m as a function of wavelength. (f,g) Performance of the metasurface for TE polarisation. 



viz. GSP-based, metasurfaces and demonstrate, both theoretically 
and experimentally, that orthogonal linear polarisations of reflected 
light can be manipulated independently, efficiently and in a broad 
wavelength range by exploiting two degrees of freedom (DOF) in the 
nanobrick geometry. 

Results 

Let us start by considering a GSP-based metasurface in a gold-glass- 
gold configuration with the top layer representing a 240-nm-period 
array of nanobricks illuminated at normal incidence (Fig. la). The 
nanobrick dimensions and Ly represent the two DOF that we are 
going to exploit, leaving all other system parameters constant (note 
that the lower gold part is considered optically thick). Considering 
the complex refiection coefficient for two orthogonal polarisations 
(denoted as TM and TE for the incident light being polarised along 
the X- or y-axis, respectively), we find immediately that, near the GSP 
resonance, the two DOF allow us to independently vary the phase of 



reflected polarisations practically within the whole phase space of 2 71 
(Fig. lb). In view of that, it is evident that full control of the reflection 
phase requires the GSP resonance to exist within the parameter space 
{L^ and Ly), otherwise, limited by the unit cell period A. The spacing 
between equidistant contour lines of the phase and the extent of 
amplitude variation are dependent on the choice of geometrical 
parameters, and can be directly related to the properties of GSPs as 
follows. Close to the GSP resonance that is characterised by strong 
absorption of the incident light (Supplementary Fig. SI), the phase 
varies rapidly over the range of n as expected for the Lorentz-like 
response of an oscillator (Fig. lb). Away from the resonance, the 
reflection phase is influenced by the retardation of the incident wave 
when being reflected primarily by the continuous gold layer (for 
smafl-area nanobricks) and by the processes of in- and out-coupling 
of the incident radiation to GSPs (for large- area nanobricks). Finally, 
for very narrow nanobricks (in the direction perpendicular to their 
propagation), the GSP gradually transforms into a surface plasmon 
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mode supported by a nanorod with the propagation constant 
increasing rapidly for smaller widths, as seen in bending of the phase 
contours and absorption locations towards shorter lengths (Fig. lb 
and Supplementary Fig. S2). 

In order to demonstrate the total control of reflected light with 
GSP-based metasurfaces, we now design three gradient metasurfaces 
that, in all cases, anomalously reflect TM waves in the same direction, 
while the reflection direction for TE waves is engineered at will. 
According to the equivalence between the generalised laws of refrac- 
tion and reflection and diffraction theory, the designed metasurfaces 
may also be viewed as flat blazed gratings, where the angle of anom- 
alous reflection is equal to the first order diffraction angle^^. 
Importantly, these blazed gratings should operate independently 
(and differently) for different polarisations, a functionality that is 
not accessible with conventional approaches. Using the grating ter- 
minology, a polarisation beam splitter (case 1) that reflects TM and 
TE waves into ± 1 diffraction order, respectively, represents a meta- 
surface with a linear phase gradient whose sign is opposite for the two 
polarisations. In this case, our goal is to implement the reflection 
coefficient r{x) = A exp{±i27ix/Asc), where A < 1 is the amplitude 
constant, A^^ is the (super cell) grating periodicity, x is the spatial 
coordinate, and the sign of the phase gradient depends on the polar- 
isation. Such a design can be derived from Fig. lb, where the 
reflection phase for both TM and TE polarisation is discretised into 
six constant-phase contour lines with a 60° step, and the appropriate 
nanobrick dimensions correspond to the intersections marked with 



circles. Note how this visual solution of the appropriate nanobricks 
contrasts to the search algorithm used in the design of birefringent 
metasurfaces in the infrared regime^^. The corresponding metasur- 
face can be constructed by sequentially placing the six elements with 
the center-to-center distance of A = 240 nm into a super cell 
(Fig. Ic), with these elements being arranged in pairs in order to 
lower the first order diffraction angle 6[^^ to 16.1°. Note that the 
described approach is based on the assumption that the coupling 
between neighbouring elements is weak, a fact that we verified by 
numerical modelling (see Supplementary Fig. S2). The cost of the 
increased super cell period is a more apparent stair-casing of the 
linear phase profile (Fig. ld,f), which results in a slightly decreased 
performance of the metasurface (Supplementary Fig. S3). Numerical 
simulations of the designed gradient metasurface show (Fig. le,g) 
that practically all reflected light for TM (TE) polarisation is con- 
tained within the + 1 ( — 1) diffraction order at the design wavelength 
of 800 nm, with the total reflectivity being limited to ~ 80% because 
of Ohmic losses in the metal. 

The selection of nanobrick dimensions for the independent phase 
control of orthogonal polarisations reUes on intersection of the cor- 
responding constant-phase lines for TM and TE polarisations, which 
is not perfect as, e.g., the 250° -phase lines do not intersect (Fig. lb). 
Ideally, one would also prefer to achieve a constant reflection ampli- 
tude, but with the two DOF in our geometry it is possible to adjust 
only two wave parameters - the reflection phases of two polarisations 
or the phase and amplitude of one polarisation^^. On the other hand. 
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Figure 2 | Functionalised gradient GSP-metasurfaces. (a-c) Metasurface polarisation beam splitter (case 2): TM (TE) waves are reflected into +1 (0) 
diffraction order, (d-f) Polarisation-independent beam steering (case 3): TM and TE waves are reflected into +1 diffraction order. (a,d) Top-view 
of designed super cells. (b,e) Theoretical performance of the metasurfaces for TM polarisation, displaying the amount of incident light reflected into the 
lowest diffraction orders m as a function of wavelength. (c,f) Performance of the metasurfaces for TE polarisation. 
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Figure 3 | Experimental verification of polarisation beam splitter (case 1): TM (TE) waves are reflected into +1 ( — 1) diffraction order. 

(a) Representative scanning electron microscopy (SEM) image of fabricated super cell, (b) Optical images of the diffraction spots at 2 = 800 nm for TM 
and TE polarisation, and reference sample (Ref.) consisting of a 100 nm gold film overlaid by 50 nm of silicon dioxide. (c,d) Measured (markers) and 
calculated (lines) diffraction efficiencies for orders Iml < 1 for (c) TM and (d) TE polarisation. 



the equivalence in operation of gradient metasurfaces and blazed 
gratings suggests that the presented approach should be robust 
against (design and fabrication) imperfections^^, an important fea- 
ture which becomes also apparent in the following consideration of 
other two cases. Keeping TM waves being diffracted into + 1 order, 
we now design metasurfaces that would reflect TE waves specularly 
(case 2) or into + 1 diffraction order (case 3). Note that the latter case 
corresponds to the polarisation-independent anomalous reflection 
(beam steering). The two different configurations correspond to 
nanobrick dimensions determined by the constant-phase line inter- 
sections in Fig. lb marked with crosses and squares, respectively. 
Clearly, the performance of case 2 for TM waves is affected by con- 
siderable variation in the reflection amplitude for different elements, 
whereas case 3 suffers from the absence of a nanobrick ensuring the 
reflection phase of 250° for both polarisations. Nevertheless, the 
designed metasurfaces in both cases perform well when being simu- 
lated as seen in Fig. 2. Note that, in case 3, instead of the pair of 
equivalent nanobricks, we used two different 250° -phase nanob ricks 
selected separately for TM and TE waves (Fig. 2d). Regarding case 2, 
it is clearly seen how TM waves are almost exclusively reflected into 
first diffraction order with a theoretical reflectivity of ~ 70% at the 
design wavelength (Fig. 2b). For TE polarisation, on the other hand, 
the incident light barely interacts with the nanobricks, which lead to a 
reflectivity resembling that of a plane gold substrate (Fig. 2c). For 
case 3, it is evident how both TM and TE waves are anomalously 
reflected in +1 diffraction order, while both polarisations dem- 
onstrate a broadband response with reflectivity of ~ 75% at >^ = 
800 nm (Fig. 2e,f). The difference between reflection into first dif- 
fraction order and the total reflectivity, which is especially apparent 
for TM polarisation (Fig. 2e), is due to reflection into higher diffrac- 
tion orders not displayed in Fig. 2. 

With the above numerical simulations illustrating the possibility 
to independently and efficiently control reflected light of orthogonal 
polarisation using gradient GSP- metasurfaces, we now move on to 
the experimental verification. Figures 3-5 display representative 
images of lithographically fabricated super cells, optical images of 
diffraction spots, and diffraction efficiencies for study case 1-3, 



respectively. In all cases, reasonable correlation is found between 
designed and fabricated super cells (Fig. 3a, 4a, 5a, and Supple- 
mentary Fig. S4), though discrepancies are also clearly visible. For 
example, in case 3, one notices how the fifth pair of elements in 
the super cell is connected to neighbouring elements along the y- 
directions (Fig. 5a). Nonetheless, the designed metasurfaces are 
robust towards fabrication imperfections, as evident from the optical 
images of the refiected light, demonstrating the independent steering 
of orthogonal polarisations (Fig. 3b, 4b, and 5b). Moreover, the per- 
formance of the GSP -metasurfaces is further quantified by the deter- 
mination of the zero and first order diffraction efficiencies (Fig. 3c,d, 
4c,d, and 5c,d). The measurements are compared with numerical 
calculations (lines) for which we use the designed dimensions of 
super cells and only one fitting parameter: the imaginary part of 
the gold permittivity is increased by four times to account for surface 
scattering and grain boundary effects in thin films as well as the 
increased damping associated with the 3-nm-thin Ti adhesion layers 
between gold-glass interfaces^ \ In general, we observe a reasonable 
agreement between the measured and calculated diffraction efficien- 
cies, verifying the independent and broadband control of TM and TE 
waves, albeit with a reduced efficiency of ~ 50% near the design 
wavelength (theory: ~ 80%, c.f. Fig. 1 and 2) related to the just 
mentioned loss channels. Measurements also show a non-vanishing 
zero order diffraction, a feature that we ascribe to the imperfections 
in the fabricated metasurfaces. 

Discussion 

We have illustrated with the examples in the previous section how 
the GSP-based metasurfaces with two DOE can be designed to inde- 
pendently, efficiently and in a broadband wavelength range control 
the phase of reflected TE and TM waves, demonstrating also the 
practical feasibility and robustness of the presented approach. 
Although the examples were limited to gradient metasurfaces with 
variation along the x-axis, it should be stressed that the reflection 
coefficients depending on two spatial coordinates can be just as easily 
implemented, thus opening up a wide range of possible functional- 
ities, including but not limited to anomalous reflection of TM and TE 
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Figure 4 | Experimental verification of independent control of TM and TE waves (case 2): TM (TE) waves are reflected into + 1 (0) diffraction order. 

(a) Comparison of designed and representative SEM image of fabricated super cell, (b) Optical images of the diffraction spots for TM and TE 
polarisation, and reference sample (Ref.) at >l = 800 nm. (c,d) Measured (markers) and calculated (lines) diffraction efficiencies for orders I ml < 1 for 
TM and TE polarisation, respectively. 



waves in two orthogonal planes, polarisation beam- splitting com- 
bined with strong focusing, and linear-to-radial (or azimuthal) 
polarisation conversion. Furthermore, we stress that the shortcom- 
ings of the nanobrick design, especially apparent in study case 3, 
can be alleviated with more complex elements with a larger num- 
ber of design parameters - one such element is the nanocross 
(Supplementary Fig. S5). Finally, we would like to emphasise that 
the manipulation of two wave parameters alternatively allow for the 



control of both the amplitude and phase of reflected light for one 
polarisation, a feature that was already used in our design of focusing 
metasurfaces to ensure the same reflection amplitude^l This idea can 
be extended further and employed to obtain a practically arbitrary 
distribution of phase and amplitude of reflected (linearly polarised) 
light. For example, one can design a metasurface acting as a beam 
splitter that splits the incident light into ±1 diffraction order by 
implementing the reflection coefficient r = A sm{2n/AscX). In this 
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Figure 5 | Experimental verification of polarisation-independent blazed grating (case 3): TM and TE waves are reflected into +1 diffraction order. 

(a) Comparison of designed and representative SEM image of fabricated super cell, (b) Optical images of the diffraction spots for TM and TE 
polarisation, and reference sample (Ref.) at 2 = 800 nm. (c,d) Measured (markers) and calculated (lines) diffraction efficiencies for orders I ml < 1 for 
TM and TE polarisation, respectively. 
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case, one should make use of the strong variation in the metasurface 
reflection ampUtude close to the GSP resonance (Fig. lb); a region 
that otherwise normally belongs to design of plasmonic super absor- 
bers^^. One should also expect that half of the incident light will have 
to be absorbed in this case (Supplementary Fig. S6). 

Overall, the presented approach establishes a new class of compact 
optical components, viz., plasmonic metasurfaces with controlled 
gradient birefringence, that has no dielectric counterparts. It can 
straightforwardly be adapted to realise new optical components with 
hitherto inaccessible functionalities as well as extended to other 
wavelength ranges. We believe that the robustness, efficiency and 
scalability of the considered configurations will stimulate the 
development of novel planar optical components compatible with 
and easily integrated in nanophotonic systems. 

Methods 

Modelling. All modelling results are performed using the commercial available finite 
element software Comsol Multiphysics, ver. 4.3a. In the calculations, it is only 
necessary to model a single unit cell or super cell by applying periodic boundary 
conditions on the vertical sides of the cell. The incident wave is in all cases assumed to 
be a plane wave propagating normal to the surface with either TE or TM polarisation 
(see Fig. la). In our gold-glass-gold configurations, the permittivity of gold is 
described by interpolated experimental values^*^ while glass, assuming to be silicon 
dioxide, takes on the constant refractive index n = 1.45. The medium above the 
nanobricks is chosen to be air. The air domain is truncated using perfectly matched 
layers to minimize reflections, while a perfect electric conductor boundary condition 
is applied on the bottom side of the optically thick gold substrate. Regarding the 
complex reflection coefficients of the GSP-based metasurfaces (Fig. lb and 
Supplementary Fig. S2), the phase is determined at the top surface of the nanobricks. 
Experimental measured diffraction efficiencies are compared with calculated values 
(Fig. 3-5) . In these calculations the imaginary part of the gold permittivity is increased 
by a factor of four to take into account additional losses related to surface scattering 
and grain boundary effects in thin films as well as the increased damping associated 
with the 3 nm Ti adhesion layers between gold-glass interfaces'^ The diffraction 
efficiencies are then obtained by calculating the reflection into the different 
diffraction orders and normalising with the reflectivity from the reference 
configuration (50 nm silicon dioxide on top of a gold substrate). 

Fabrication. The GSP-based metasurfaces are fabricated using electron-beam 
lithography and lift-off technique. On a base substrate of silicon, 3 nm of titanium, 
100 nm of gold and 3 nm of titanium are sequentially deposited by E-beam 
evaporation and subsequently 50 nm of Si02 is deposited by RF- sputtering in a 
Cryofox 600 system equipped with a Temescal SuperSource2 E-beam source and a 
MAK magnetron RF-sputter source. Following a spin coating of 120 nm PMMA 
950K A2, the gold nanobrick structures with an interparticle distance of 240 nm are 
lithographically defined by E-beam lithography and formed by E-beam evaporation 
of 3 nm Ti and 50 nm gold followed by a lift-off process. One should note that the 
three layers of 3 nm titanium are used as adhesion layers at the silicon-gold and glass- 
gold interfaces. The nanobricks are arranged in 28.8 fim X 28.8 fim arrays comprised 
of 10 and 100 super cells in the x- and y-direction, respectively, and they are exposed 
at 30 kV using a JEOL-640LV electron microscope equipped with an ELPHY 
Quantum lithography system. The resulting structures are imaged using a Hitachi S- 
4800 field emission scanning electron microscope. 

Optical characterization. Experimental characterisation of the fabricated structures 
is performed by directly measuring the intensity of the anomalously reflected laser 
beam into various diffraction orders (0 and ± 1). A low-power (~5 mW) continuous- 
wave laser beam from Spectra- Physics 3900 S Ti:Sapphire laser, tunable in the 
wavelength range 700-860 nm, was weakly focused using a 5X objective into a spot 
with diameter ~ 20 fim, so that it was illuminating an array of gold nanobricks, not 
extending beyond them. A beam splitter right before the surface of the sample allowed 
separating the incident and the reflected signals and measuring the intensity of the 
latter, which is then normalised to the intensity of the beam reflected from a smooth 
surface without nanobricks. In order to visualise the effect, we also projected the 
reflected pattern onto a charge -coupled device camera (Fig. 3b, 4b, and 5b). 
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